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SUMMARY 


Al 

The  flight-test  techniques  for  measuring  continuous  atmospheric 
turbulence  are  briefly  outlined  and  available  measurements  of  the  power 
spectrum  of  atmospheric  turbulence  reviewed.  The  measured  results  are 
examined  in  order  to  determine  how  well  they  approximate  the  idealized 
models  of  a  stationary  Gaussian  random  process  and  of  isotropic  turbu¬ 
lence.  Finally,  some  recent  results  on  the  variation  with  altitude  of 
the  probability  distributijn  of  the  root -mean -square  gust  velocity  arp 
presented  and  the  application  of  these  results  to  operational  response 
history  calculations  indicated. 
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NOTATION 


an  acceleration 
*  see  Equation  (14)  in  imin  text 

a,  airplane  longitudinal  acceleration 

c  average  chord,  ft 

KK.s)  airplane  gust  response  function,  Equation  (15)  in  main  text 

K  4W 

gnpS£ 

h(t)  unit  inpulse  response 

H(oj)  frequency  response  function 

1  distance  from  vane  to  airplane  center  of  gravity,  ft 

L  scale  of  turbulence,  ft 

■  slope  of  lift  curve  per  radian 

M(a„)  average  number  cl  maximum*  per  second  exceeding  given  value  of 
ajj  in  operations 


N(  ) 

°i.°t 

P(cr) 

$«v 


s 


t 

u 

v 

V 

w 


characteristic  frequency,  Equation  (11a)  in  main  text 


average  number  of  maximums  per  second  exceeding  given  va) 
specified  argument  for  Gaussian  disturbances 


ue  of 


output  responses 
probability  distribution  of  a 
probability  distribution  of  ctw 
wing  area,  sq  ft 
time,  sec 

longitudinal  gust  velocity,  ft/sec 
lateral  or  side  gust  velocity,  ft/sec 


airplane  forward  speed,  ft/sec 
vertical  gust  velocity,  ft/sec 


airplane  weight,  lb 


ft 


»a  airplane  normal  or  vertical  velocity,  ft /sec 

a  angle  of  attack,  radians 

y  flight-path  angle,  wa/V 

6  pitch  attitude,  radians 

6  pitch  velocity,  radians/sec 

A.  gust  wave  length,  ft 

p  air  density, slugs/cu  ft 

cr  root-mean-square  deviation 

$(nj)  power-spectral-density  function 

H  reduced  frequency  (oj/V),  radians/ft 

cc  frequency,  radians/sec 


Subscripts 
i  input 

o  output 
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ATMOSPHERIC  TURBULENCE  ENVIRONMENT  WITH  SPECIAL 
REFERENCE  TO  CONTINUOUS  TURBULENCE 

Harry  Press* 


*  1 .  INTRODUCTION 

Aeronautical  studies  of  atmospheric  turbulence  have,  in  the  past,  been  largely 
based  on  simplified  concepts  involving  discrete  gusts  of  idealized  shapes  such  as 
sharp-edged  gusts,  triangular  gusts,  and  ramp-platform  shaped  gust  profiles.  During 
the  last  few  years,  a  number  of  developments  have  contributed  to  making  a  more 
detailed  and  realistic  approach  to  the  rough-air  problem  more  urgent.  These  include 
on  the  one  hand,  the  increased  importance  of  airplane  response  to  gusts  in  stability’ 
ant  guidance  problems.  Or^  the  other  hand,  the  need  for  the  rapid  exploitation  of 
perforBance  capabilities  has  served  to  place  greater  emphasis  on  a  more 
efficient  structural  analysis.  As  a  consequence,  the  gust  response,  which  was  formerly 
of  serious  concern  only  for  large  transport -type  airplanes  and  bombers,  has  also 
become  of  concern  for  a  much  wioer  range  of  airplane  types  as  well  as  for  missiles, 
or  una  e  y,  during  this  period,  new  mathematical  and  experimental  tools  have  been 
developed  which  served  to  make  a  more  detailed  study  of  atmospheric  turbulence 
feasible.  Noteworthy  developments  occurred  in  the  theory  of  random  processes,  in  the 
statisticai  theories  of  turbulence,  and  in  methods  for  the  measurements  of  turbulence 

in  flight,  and  these  developments  have  served  to  provide  the  basic  tools  for  a  deeper 
approach  to  gust  problems.  aeeper 

As  a  consequence  of  these  new  developments,  the  last  5  years  have  seen  a  rapid 
dwoiopmont  in  the  study  of  continuous  atmospheric  turbulence  and  airplane  benavior 
in  rough  air.  The  detailed  characteristics  of  continuous  atmospheric  turbulehce 
have  been  measured  in  a  number  of  studies.  These  measurements  have  served  to  pro¬ 
vide  a  new  insight  into  airplane  behavior  under  continuous  rough-air  conditions. 
References  1  -  25  list  some  of  the  more  important  contributions  in  these  areas. 

is  °f  detailed  “^urements  of  continuous  atmospheric  turbulence 

s  still  quite  limited,  sufficient  data  appear  to  be  available  to  warrant  a  review 

lst^s8o^tZZr?c  r*  le&8t  &  Pr8lifflinary  ^s^ment  of  the  general  character¬ 
istics  of  atmospheric  turbulence.  In  the  present  paper,  the  general  methods  used  In 

e  measurement  of  continuous  atmospheric  turbulence  in  flight  will  be  reviewed 

information  “  n'hf  8mmirlze  80ms  of  tbe  more  significant  features  of  the 

information  available  on  the  characteristics  of  atmospheric  turbulence.  In  order 

to  provide  a  compact  description  of  turbulence  measurements  and  to  provide  a  guide 

to  the  selection  of  simple  models  for  analytic  studies,  some  of  the  basic  concepts 

and  results  of  random  process  theory  and  the  statistical  theories  of  turbulence  will 

be  described  in  the  appendix  to  this  paper.  The  important  concepts  of  a  staUonlry 

fh  “8;  “?,*  GaU88ian  random  process  w111  »>«  defined,  and  the  characteristics 
f  isotropic  turbulence  will  be  described.  These  Idealized  models  are  used  in  the 

assessment  of  the  experimental  measurements  of  turbulence  in  the  main  body  of  this 
paper- 
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In  addition  to  the  foregoing  material  which  is  confined  to  the  local  characteristics 
of  atmospheric  turbulence,  consideration  will  also  be  given  to  some  of  the  broader 
climatological  variations  in  turbulence.  In  particular,  some  recently  obtained 
estimates  of  the  probability  distribution  of  the  root-mean-square  gust  velocity  and 
its  variation  with  altitude  are  presented  and  discussed.  The  application  of  these 
results  to  the  problem  of  estimating  gust  and  gust  load  histories  in  airplane  opera¬ 
tions  is  also  considered  briefly. 


2.  METHODS  OF  MEASURING  ATMOSPHERIC  TURBULENCE  IN  FLIGHT 

The  need  for  the  consideration  of  the  continuous  character  of  turbulence  has,  in 
recent  years,  led  to  the  development  of  new  and  improved  methods  of  measuring 
turbulence  in  flight.  During  the  last  few  years,  two  significant  contributions 
were  made  in  this  area.  These  are:  first,  the  general  approach  pioneered  by 
Clementson,  at  Massachusetts  Institute  of  Technology1  for  the  determination  of 
power  spectra  of  atmospheric  turbulence  by  means  of  power  spectral  techniques;  and 
secondly,  the  subsequent  development,  at  both  the  Cornell  Aeronautical  Laboratory21* 24 
and  the  N.A.C.A.  ll«  23  of  direct  turbulence  measurements  by  means  of  flow  direction 
vanes.  These  techniques  have  been  used  to  provide  basic  data  on  the  detailed 
characteristics  of  turbulence.  A  few  remarks  on  these  methods  and  their  limitatipns 
thus  appear  appropriate  in  order  to  appreciate  the  limitations  of  the  available 
information^ 


2.1  Spectral  Technique 

In  Reference  1,  Clementson  first  applied  the  concepts  of  random  process  theory 
(see  Appendix)  to  the  problem  of  airplane  behavior  in  rough  air.  *  In  this  approach, 
the  power  spectra  of  the  random  disturbance  and  of  the  airplane  response  play  a 
central  role.  Clementson  showed  how  output  response  measurements  of  a  linear  air¬ 
plane  system  could  be  used  to  determine  the  power  spectrum  of  atmospheric  turbulence. 
The  approach  used  is  based  on  the  application  of  the  relation  between  the  power 
spectrum  of  a  disturbance  and  the  power  spectrum  of  the  response  to  the  disturbance 
for  a  linear  system.  The  relation  between  the  i>ower  spectra  is  given  by 

$oe<n)  =  Ssi(oj)  j  H  ('>))  I  2  (1) 


where 

w  =  circular  frequency,  radians/sec  J 

$0(<xi)  =  the  power  spectrum  of  the  response  or  output  0(t) 

=  the  power  spectrum  of  the  input  or  gust  disturbance  I(t) 

H(oj)  =  the  frequency  response  function  for  the  system  and  in  the  present  case 

defines  the  response  of  the  system  to  unit  sinusoidal  gusts  at  the  various 
frequencies. 

The  application  of  this  relation  to  the  measurement  of  gust  spectra  is  schematically 
illustrated  in  Figure  1.  The  reduced  frequency  argument  fl  (  =oj/v)  which  is  equal  to 


™rlations  oYL":  Tf  'lvele"8th'  18  used  for  «»  abscissa  Inasmuch  as  the  space  3 
iations  of  the  turbulence  are  of  principal  concern  in  aeronautical  applications 

0llne'>ZmirlccelerTntS  ‘In  P0'"ir  SPe°tr”'  °f  “  ,Irpl“e  response  8'Kh  18  «ir-‘ 
teristlcs  anffhe  a  „  Se00nd  d“1've  repre8“ts  the  airplane  response  charac- 
terlstics.  and  the  sketch  shorn  In  typical  of  the  nornal  acceleration  response  at 

the  center  of  gravity  of  some  current  airplane  types.  The  frequency  response 
function  normally  contains  several  peaks  .hich  reflect  the  varloos  airplane  response 

period  JZ  »M  h  T'  !  rSt  Peak  18  ‘  reflecti0"  of  the  airplane  short- 
period  mode  which  normally  peaks  at  around  1/2  cycle/sec,  or  a  corresponding  gust 

Xett„  it  2:  t:  i1'?00  ft- The  second  8hd*-  -  the3ch  ise r 

fto  5  cvolf/  ?  ln  Sy*~  rlcal  bendin«  which  frequently  is  present  at 

5  cycles/sec  and  corresponds  to  gust  wavelengths  of  the  order  of  100  ft  Thus 

wavelengths  of  from  less  than  100  ft  to  several  thousand  feet  are  oJ  principal 

concern  in  many  aeronautical  applications.  Special  problems,  such  as  missile 

^i«tpcr:n;  TtTVi  ‘TJr  °onsider8ti°" « i°n«er 

souared^fThe  i  ?  '  by  dl’wln*  the  response  spectre,  by  the  «*>lltude 

r,6Bo"se  ,-ti-  ^  «» »«•«■  wt. 

min1!  I*  Clear  fr0!"  Equatlon  (I)  that  the  P^er  spectrum  of  the  input  may  be  deter¬ 
mined  from  measurements  of  the  power  spectrum  of  an  airplane  response  JlLe  if  the 

frequency  response  function  for  the  airplane  is  known.  This  is  essentially  the 
procedure  first  used  by  Clementson  and  subsequently  hy  others.  In  these  applications 
he  general  practice  has  been  to  determine  the  frequency  response  function  by 

SS  °f  C0UrSC'  ‘1S0  ^  U8ad  t0  calculate  ^ponses  ^ 

and  has  found  Tie  lagT*  **"  C@ntral  006  aer0MUtic*1  applications 

SJT  ‘"™8tl«atora  "“‘ably  Lappi,  at  the  Cornell  Aeronautical  laboratory,  and 
““““"“'‘t8  iMtltnte  of  Technolonr’.  hare  extsnded  this  approach  to 

tte  ™tl«l’^  hT“  df 'r“nln‘  “T'r*1  hhXWhht.  of  the  turbulence  such  as 

the  vertical  and  head-on  or  longitudinal  components.  This  extens<on  requires  the 

affect^Tth^^f  °fh8®Veral  0Utput  resPonse*.  one  being  significantly 
airected  by  the  several  turbulence  components. 

JZ  ,t”d‘"t“rd“Ce?  °r  T  CCW,"Cnt!'  lL(t)  *nd  «*  relatione  for  the 

uukyuc  responses  are  given  by 


(2) 


V'1  =  -  t.idt,  +  rhu(t1ni,t .  tiMti 

°!<t>  =  Oil**!*1!*1  -  ‘iM‘t  t  Ovwt  -  tjldt, 

ahere  h,,(t>  la  the  unit  inpulse  reaponse  In  Oj  to  the  disturbance  I,.  The  relations 
?!? ?_  “d  cr08s  hi*81™  <•«  Appendix)  of  the  disturbance.  and 
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where 


the  power  spectra  of  the  input  disturbances  if  i 
between  the  input  disturbances  if  i  ^  j 


-  j  or  the  cross  spectra 


the  power  spectra  of  the  outputs  if  i  =  j  or  the 
the  outputs  if  i  £  j 


cross  spectra  between 

* 


HijM  =  frequency  response  function  defining  the  response  in  0.  to  unit  sinusoidal 
disturbances^ in  Ij  at  the  various  frequencies  and  is  given  by 

Hij(<y)  =  ighjj  (t)e"^dt.. 


The  asterisk  designates  the  complex  conjugate.  If  the  Hh*s  are  known 

irn  1  iiAn  a  J?  1L.  <r\  •  II  W 


values  of  the  I  _  __  __ 

spectra  and  cross  spectra 
forward.  J 


measured 


lj  s  may  be  used  in  Equation  (3)  to  determine  the  required  disturbance 
T!le  extension  to  three  or  more  components  is  straight- 


The  principal  difficulty  with  this  approach  for  the  determination  of  atmospheric 
turbulence  spectra  is  the  requirement  for  reliable  estimates  of  the  airplane  fre¬ 
quency  response  functions  over  the  wide  range  of  frequencies  which  are  normally  of 
f*  ffquency  Unctions  can  generally  be  obtained,  but 

aooroJh  Jo  ^  ,6t  Considered  •de«U6t«  to  warrant  the  general  use  of  this 

approach  for  the  present  purpose.  Experimental  methods  of  determining  the  frequency 

response  functions  of  an  airplane  for  the  gust  condition  appear  to  be  promising 

is  Tt  r*  Z  lr\l  ^  a0  etrly  8tat®  °f  develoixoent*  An  additional  complication 
is  introduced  by  the  unavoidable  presence  in  flight  tests  of  pilot  control  motions 

of  thTairnr"06!-  ^  dir6Ct  °f  "ea8urlng  turbulence  that  are  independent’ 

of  the  airplane  frequency  response  functions^  desirable.  The  flow  direction  vane 
echnique  described  in  the  next  section  appears  to  meet  this  need. 


2.2 


Flow  Direction  Vane  Technique 


.  T!*e  U®e.°!  flow  direction  vane  measurements  for  the  study  of  turbulence  was 
CrMe°H?  t^TcT“.S  C°^e11  Aer0naUtical  tAboratory21’  24  and  by  Chilton  and 

the  flow  direrHnn  '  gUre  2  illustrates  schematically  the  relations  between 

the  flow  direction  and  the  gust  velocity  for  an  airplane. 

The  airplane  angle  of  attack  a  is  given  by 


a. 


=  9  +  y 


where 


(4) 


0  ~  the  pitch  attitude,  radians 

y  ~  wa/v 

wa  =  the  airplane  normal  velocity,  ft/sec 

V  =•  the  airplane  forward  speed,  ft/sec. 


ni  i  WMir«lfatmn4MWii*r  1  . .  <•*•**■  ■*** 


5 

* 

For  small  disturbances,  the  angle  sensed  by  the  vane  ay,  neglecting  the  upwash,  is 
then  given  by 


°y 


+  aD 


(5) 


where 


wg  =  the  normal  (or  vertical)  gust  velocity,  ft/sec 
l  =  the  distance  from  the  vane  to  the  airplane  center  of  gravity. 

* 

From  Equations  4  and  5,  the  vertical  gust  velocity  is  then  given  by 

wg  =  VfOy  -  6)  -  wR  +  id  (6) 

Equation  6  indicates  that  the  time  history  of  gu3t  velocity  sensed  by  the  airplane 
can  be  determined  from  the  simple  addition  of  measurements  of  airplane  vane  angle, 
pitch  attitude,  airplane  vertical  velocity,  end  pitch  velocity.  The  use  of  the 
equivalent  expression  for  the  reduction  of  sideslip  angle  measurements  to  obtain 
the  side  gust  component  does  not  appear  to  offer  any  additional  serious  difficulties. 

In  principle,  the  application  of  Equation  (6)  to  gust  measurements  is  straightforward. 
In  practice,  a  number  of  difficulties  are  encountered  which  tend  to  limit  the  accuracy 
and  the  range  of  frequencies  which  can  be  covered.  These  difficulties  stem  essentially 
from  limitations  in  measurement  techniques.  As  might  be  expected,  the  instrument 
accuracy  requirements  for  such  applications  are  quite  severe.  For  an  airplane 
flying  at  300  ft/sec,  a  1-ft/soc  vertical  gust  gives  rise  to  an  angle-of-attack 
change  of  1/300  of  a  radian,  or  roughly,  0.2°.  It  is  thus  clear  that  in  order  to 
achieve  an  accuracy  of  1  ft/sec,  it  is  necessary  to  measure  the  required  angles  to  al 
least  an  accuracy  of  0.1°  or  0.05°.  This  high  degree  of  sensitivity  is  difficult 
to  achieve  but  can  be  approximated  in  practice  when  some  precautions  are  taken. 

In  addition,  good  frequency  response  characteristics  to  relatively  high  frequencies 
are  required  in  order  to  avoid  distortions  due  to  magnifications  and  phase  shifts. 


Flow  direction  vanes  with  good  frequency  response  characteristics  do  not  appear 
difficult  to  obtain.  Standard  metal  vanes  frequently  have  response  characteristics 
that  are  essentially  flat  to  10  cycles/sec.  Lightweight  balsa-wood  vanes  with  good 
response  to  much  higher  frequencies  are  not  difficult  to  obtain.  Angle-of-attack 
■  pressure  sensors  have  also  been  considered.  Perhaps  the  most  troublesome  of  the 

measurements  is  the  measurement  of  pitch  attitude.  Although  good  frequency  response 
characteristics  can  be  obtained,  conventional  pitch  attitude  gyro  recorders  are 
subject  to  low  frequency  drift.  ThiB  difficulty  can  largely  be  avoided  by  obtaining 
pitch  attitude  from  the  integration  of  pitch  velocity  measurements.  Some  investiga¬ 
tors  have  had  recourse  to  the  use  of  a  sun-camera  technique  for  studies  covering 
very  low  frequencies23.  Airplane  vertical  velocity  is  generally  difficult  to  measure 
directly  but  the  integration  of  normal  acceleration  measurements  appears  to  provide 
a  satisfactory  measure  of  vortical  velocity  when  the  initial  conditions  are  adjusted 
,  through  the  use  of  a  sensitive  pressure  altimeter  or  statoscope. 


SSKVJPbA 
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An  additional  practical  difficulty  is  worth  noting.  The  necessity  for  getting 
the  vane  [Measurements  out  of  the  airplane  flow  field,  in  general,  requires  the  use 
of  a  vane  mounted  on  a  boom.  Booms  mounted  at  the  airplane  center  line  have 
appeared  best  for  rigidity  reasons  and  have  the  additional  advantage  of  avoiding 
the  complications  of  the  effects  of  rolling  and  yawing  motions.  Boom  lengths  of  the 
order  of  one  chord  length  have  been  used  and  appear  to  be  satisfactory.  However, 
in  practice,  the  boom  stiffness  is  generally  not  as  good  as  is  desirable.  Boom 
vibrational  frequencies  thus  become  e  problem,  but  it  is  generally  possible  to  keep 
the  boom  frequencies  sufficiently  high  (5-15  cycles/sec)  to  be  out  of  range  of  the 
more  significant  gust  frequencies  which  are  generally  at  the  lower  frequencies. 

The  sensitivities  used  in  recording  the  various  quantities  can  considerably 
influence  the  reliability  of  the  spectrum  measurements.  Because  of  the  general 
character  of  the  gust  spectrum  and  the  spectra  of  the  various  other  airplane 
measurements,  the  measurements  will,  in  general,  contain  yepy  little  power  at  the 
higher  frequencies.  Thus,  the  results  obtained  at  the  higher  frequencies  particu¬ 
larly  are  susceptible  to  large  errors.  As  a  consequence,  studies  aimed  at  covering 
a  wide  range  of  frequencies  may  require  different  sensitivities  and  particularly 
higher  sensitivities  for  the  higher  frequency  ranges.  The  fact  that  w  and  6 
are  determined  from  the  integration  of  acceleration  and  pitch  velocity  measurements 
makes  the  low  frequency  accuracy  of  w&  and  Q  critical  since  small  errors  in  say  6 
(or An) at  low  frequencies  yield  larger  errors  in  derived  values  of  9  (or  w  ). 

A 

The  numerical  determination  of  the  pcwer  spectra  from  the  time  history  data  is, 
in  general,  a  task  of  some  delicacy.  Care  is  required  in  the  choice  of  reriing 
intervals,  and  length  of  sample  record  in  order  to  obtain  reliable  results.  The 
actual  numerical  operations  in  the  determination  of  the  auto-correlation  function 
and  power  spectrum,  as  defined  in  Equations  (3)  and  (4)  of  the  Appendix,  must 
also  be  carefully  planned  to  avoid  distortions  arising  from  the  limited  data  and 
from  the  numerical  operations  themselves.  These  problems  are  particularly  important 
in  the  gust  spectra  determinations  and  are  discussed  in  some  detail  in  Reference  20. 

The  full  expression  given  by  Equation  (6)  is  not  required  over  the  full  frequency 
range  of  interest  which  covers  the  wavelength  range  of  1  ft  to  perhaps  50,000  ft. 

In  portions  of  this  wide  range  of  wavelengths  or  frequencies,  some  of  the  terms 
become  negligible.  At  high  frequencies,  say  above  2  or  3  cycles/sec  for  many 
current  airplane  types,  the  airplane  pitch  and  vertical  motion  response  are  generally 
small  relative  to  the  gust  angle  and  can  be  neglected.  For  this  case,  Equation  (6) 
reduces  to 


»g  =  VXy  (7) 

and  the  vane  measurement  alone  provides  a  direct  measure  of  the  turbulence^  At  very 
low  frequencies,  generally  less  than  about  0.25  cycles/sec,  the  airplane's  inherent 
stability  acts  to  keep  the  airplane  at  a  fixed  angle  of  attack,  and  the  vane 
measurements  are  then  of  r.o  use.  The  gust  velocity  for  these  low  frequencies  is 
then  given  by 

**  =  -v9  ■ ». 
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and,  in  general,  depends  upon  small  differences  between  the  vertical  motion  and  the 
pitching  motion.  The  accuracy  requirements  on  these  two  measurements  consequently 
become  more  critical.  The  reliability  of  conventional  pitch-rate  gyros  is  somewhat 
questionable  at  these  lower  frequencies,  and  pitch  attitude  measurements,  using  a 
sun-camera  technique,  have  for  this  reason  been  used  in  this  frequency  region23. 

In  order  to  provide  an  indication  of  the  magnitude  of  the  corrections  introduced 
by  the  various  airplane  motions  in  the  airplane  short-period  frequency  region,  some 
recent  results  obtained  in  connection  with  a  study  of  the  behavior  of  a  large 
swept-wing  airplane  in  rough  air  are  summarized  in  Figure  3.  Figure  3  shoro  the  basic 
spectrum  of  \fay  (the  vane  indicated  gust  velocity).  The  modified  gust  velocity 
spectra  obtained  when  the  corrections  for  vertical  motions  ( Va^  -  w&)  and  when  the 
corrections  for  vertical  motions  and  pitching  motions  (VOy  -  vd-  wa  +  1$)  are 
included,  are  also  shown. 

Figure  3  indicates  that  the  corrections  for  the  airplane  motions  tend  to  increase 
the  gust  spectrum  at  very  low  frequencies  (Cl  <0.003)  indicating  that  the  airplane 
pitching  and  vertical  motions  act  to  reduce  the  airplane  angle  of  attack  due  to  the  *- 
gust.  At  the  higher  frequencies  (0.003  <  Cl  <  0.010)  the  airplane  motions  increase 
the  airplane  angle  of  attack  above  the  gust  angle-of-attack  change.  This  is  largely 
a  consequence  of  the  pitching  motions  in  the  short-period  frequency  region.  The 
magnitude  of  the  corrections  is  quite  large  in  this  region  reducing  the  gust 
spectrum  to  about  one-third  that  indicated  by  the  vane  angle  alone.  At  still  higher 
frequencies,  the  difference  between  the  vane  indicated  gust  velocity  and  the  full 
expression  is  seen  to  be  small,  although,  in  the  present  case,  the  situation  is 
complicated  by  the  presence  of  vibration  effects  associated  with  the  aiVpl&ne 
fundamental  vibration  mode  which  has  a  frequency  corresponding  to  about  O  =  0.012 
radians  per  ft.  It  is  thus  clear  from  the  foregoing  that  the  correction  for  pitching 
and  vertical  motions  is  quite  large  in  the  airplane  short-period  region  and  is 
required  for  reliable  gust  spectra. 

In  spite  of  the  many  difficulties  encountered  in  practice,  the  flow  direction 
technique  for  turbulence  measurement  appears  to  provide  the  most  reliable  method 
so  far  devised  for  the  measurement  of  atmospheric  turbulence  in  flight.  The 
experimental  results  to  be  described  subsequently  will  be  largely  restricted  to 
such  measurements. 

2.3  Other  Methods 

Two  other  methods  of  measuring  turbulence  deserve  mention  and  will  be  discussed 
briefly. 


Inasmuch  as  pressure  measurements  in  flight  are  a  well-developed  art,  the  high- 
frequency  fluctuations  of  a  sensitive  and  high-frequency  airspeed  system  provide  a 
useful  method  of  measuring  the  longitudinal  (head-on)  component  of  turbulence.  The 
longitudinal  gust  velocity  u(t)  would  appear  to  be  given  simply  by 

u(t)  =  A V(t)  -  r  ax(t1)dt1  (9) 
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where 

AV  =  the  fluctuation  of  the  measured  airspeed,  ft/sec 
ax  =  ^*e  l°n8itudinal  airplane  acceleration,  ft/sec. 


Because  of  the  high  airplane  inertia,  the  acceleration  ax  is  frequently  negligible 
at  the  higher  frequencies  and  the  fluctuations  of  the  airspeed  may  sometimes  be 
taken  as  a  direct  measure  of  the  turbulence.  For  many  airplanes,  this  appears  to  be 
reliable  for  gust  frequencies  down  to  at  least  1  cycle/sec.  The  longitudinal  gust 
velocities  have  not.  however,  yet  received  much  detailed  consideration  inasmuch  as 
their  effects  on  the  airplane  are  generally  small. 

A  large  body  of  atmospheric  turbulence  measurements  has  been  collected  in  recent 
years  by  meteorologists  using  wind  measuring  equipment  mounted  on  towers  These 
measurements  are.  of  course,  restricted  to  the  lower  several  hundred  feet  of  the 
atmosphere.  Also,  the  direct  applicability  of  these  measurements  to  aeronautical 
questions  appears,  at  this  time,  to  be  open  to  some  question,  and  for  these  reasons 
as  well  as  .or  space  limitations,  it  was  decided  to  make  no  direct  use  of  this 

material  in  the  present  study.  Some  qualitative  references  to  these  results  xiU 
however,  be  made.  -  ' 


3.  CHARACTERISTICS  OF  ATMOSPHERIC  TURBULENCE 
3.1  General  Remarks 


In  this  section  of  the  paper,  some  of  the  significant  features  of  the  available 
information  on  atmospheric  turbulence  are  reviewed.  No  effort  will  be  made  to  be 
exhaustive  ih  covering  the  available  data,  but  rather  the  intent  will  be  to  present 
illustrative  results  cn  the  various  questions  of  interest.  The  discussion  will  be 
divided  into  two  parts,  one  concerned  with  what  might  be  called  the  local  detailed 
characteristics  of  atmospheric  turbulence,  the  other  concerned  with  the  broader 
climatological  variations  of  turbulence  intensity. 


In  considering  the  detailed  statistical  characteristics  of  atmospheric  turbulence, 
it  will  be  helpful  to  determine  how  well  the  measurements  can  be  approxinated  by 
simple  theoretical  models  that  have  been  developed  both  in  random  process  theory 
and  in  general  theories  of  turbulence.  Such  evaluations  are  important  both  for 
guiding  analytic  studies  of  airplane  behavior  which  can  be  simplified  hy  the  use  of 
idealized  models  and  also  for  planning  turbulence  measurement  surveys.  The  specific 
turbulence  properties  that  are  pertinent  and  are  of  present  interest  are: 

(a)  The  power  spectra  of  the  component  velocity  fluctuations 


(b)  Homogeneity  of  the  turbulence 


(c)  Stationarity  of  the  turbulence 


(d)  Gaussianness  of  the  velocity  fluctuations 

(e)  Isotropy. 
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These  properties  are  defined  analytically  m 
6  ^  "  “  brl6fly  ^  Sl-1—  and6 

provides  ^detailed 'descriptiOT^fhthe^ta  °^s  t  i  ra?d°f  disturbance  in  general 

and.  in  particular,  describes  the  frequency ‘contenV Char8cteristlcs  o{  the  disturbance, 
case  °f  velocity  fluctuations  in  turbulence  with  whVh  dlstur*>ance.  For  the 
spectrum  describes  how  the  kinetic  energy  is  distrfh^  H*  ***  concerned«  the  power 
wavelengths.  Inasmuch  as  airplanes  •»  huted  with  frequency  or  gust 

resDon l0US  *avelength8'  this  representations*  wrt ic Jla  * degress  to  turbulence 
response  studies.  ls  Particularly  useful  in  aeronautical 

In  the  statistical  c^terls^cr^fX^SJen1"1’^  h  Peralstenc*  “<1  !n*arl«nee 
respect!™,,  Thus,  lor  .  realon  o?  th  -  time  displacements, 

* *  airplane  measurements  would  yield  the  same  . “rb“lence*  Jt  "<>uld  be  expected  that 
starting  point  of  the  measurement,  ^^istics  regardless  of  L 

be  invariant  for  different  times  of  the  su^l  ^  **  ““™«t.  would 

involving  stationary  and  homogeneous  disturbance  T  8Pecial  case  of  Problem 
and  has  received  considerable  attention.  Particularly  amenable  to  analysis 

associated  ,1th  the  dlstnrbaoce.  This  s£ciflc«.  t  ‘>UltT  dl»trlbutlons 
derivations  In  application  of  specific  1”  tremendous!,  slwllfies  the 
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characteristics  of  the  turbulence  measurad  „!?  a*?11  that  the  st»tiatical 
rection.  A  further  important  consequence  of  iLtm?!  °f  ^  *irpIane  «i*ht 
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3.2.1 


Local  Turbulence  Characteristics 
Time  History  Characteristics 


atmospheric  turbull^e  characteristics  of 

concurrent  time  histories  of  the  vert  Wi  ^  VV  lghtf  Pigure  4  Presents  short 
J-»  For  comparlso  the  ^l.t  ^Ur  *"d  **  sW«  07  lateral 

la  S0»tlme.  used  as  .  measure  of  thH^w!  ?  ',0r",I  ““‘•"‘lo.,  .hid. 

The  first  point  of  Interest  ls  the  Irregularis’  '  s*°*"  ln  ttK  figure, 

fluctuations,  resembling  a  random  type  of  oscillat!™10  °har*cter  of  the  turbulence 
■a  dhite  tenet..  Is  of  oonrss  the  basic  re^^he 
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tion  of  statistical  techniques.  Although  there  is  some  tendency  for  the  normal 
acceleration  to  follow  at  least  the  more  rapid  vertical  gust  velocity  changes,  the 
two  records  are  very  different.  Perhaps  the  major  difference  noticeable  to  the  eye 
is  the  absence  of  the  long  waves  in  the  acceleration  record.  This  difference  results 
from  the  airplane's  lack  of  sensitivity  to  these-  longer  gust  wavelengths.  Another 
point  of  interest  is  the  lack  of  correlation  between  the  vertical  and  side  gust 

velocities. 


3.2.2  Power  Spectrum  of  Turbulence 

A  number  of  flight  measurements  of  the  spectrum  of  atmospheric  turbulence  have  now 
been  obtained  at  Corne It  Aeronautical  Laboratory,  Massachusetts  Institute  of  Tech¬ 
nology.  and  the  N.A.C.A.  Additional  studies  are  currently  under  way.  These  measure¬ 
ments  reveal  a  rather  surprising  degree  of  uniformity  in  the  general  shape  of  the 
spectrum  For  this  reason,  no  effort  will  be  made  to  present  the  many  measurements 
available  but  rather  one  set  of  measurements  for  vertical  gust  velocity  is  shown 
in  Figure’ 5.  This  measured  spectrum  was  obtained  in  Reference  23  and  is  essentially 
typical  of  most  of  the  measurements  available  and  has  the  added  distinction  of 
covering  at  one  time,  a  very<»wide  range  of  gust  wavelengths,  the  range  of  gust 
wavelengths  covered  being  from  10  ft  to  about  50.000  ft.  The  result  shown  in  the 
figure  for  the  spectrum  is  given  in  three  overlapping  sections  corresponding  to  the 
reduction  procedure  used  for  the  several  frequency  ranges  as  given  in  Equations  (0) 
to  (8).  The  most  significant  aspect  of  the  spectrum  of  the  turbulence  is  the  rapid 
decrease  with  increasing  frequency  of  the  spectral  power.  In  fact,  for  frequencies 
greater  than  D  =  0.001,  the  frequency  decreases  approximately  as  w  ,  or  perhaps 
within  the  accuracy  of  the  data  as  (Ts/3.  a  rate  predicted  by  the  theory  of  isotropic 
turbulence  At  the  lower  frequencies,  fi  <  0.001,  a  tendency  is  indicated  in  these 
measurements,  as  in  other  studies,  of  a  flattening  of  the  spectrum.  This  region 
of  frequencies  has.  however,  so  far  received  little  attention. 

In  many7  analytical  studies,  the  following  analytic  representation  for  the  spectrum 
of  vertical  (or  lateral)  gust  velocity  has  been  used.  The  expression  is  given  by 


$(H)  =  o', 


2  L_  1  ±  3ffL2 


*  77(1  +  q2l2)2 


(10) 


where 

crw  =  the  root -me an- square  gust  velocity,  ft/sec 
L  =  the  scale  of  turbulence,  ft. 


The  quantity  cr  provides  a  simple  measure  of  the  magnitude  of  the  spectrum  or  the 
intensity  of  the  turbulence,  while  the  scale  of  turbulence  L  provides  a  measure  of 
the  average  eddy  size  and  indicates  the  frequency  at  which  the  spectrum  changes 
from  a  flat  spectrum  to  one  which  decreases  rapidly.  This  frequency  is  given 
approximately  by  0  =  1/L.  The  curve  for  <rw2  =  33  and  L  =  1,000  ft  is  shown  in  Figure  5 
for  comparison  with  the  measurements  and  is  seen  to  yield  a  good  approximation  to 
the  measurements  over  almost  the  entire  range  of  frequencies  covered. 
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Except  for  differences  in  intensity,  or  in  crm2,  the  result  shown  in  Figure  5  for 
the  spectrum  of  vertical  gust  velocity  appears  to  be  representative  of  the  measure¬ 
ments  obtained  in  flight  tests.  However,  it  should  be  noted  that  the  frequency 
ranges  covered  in  most  studies  are  quite  limited.  Also,  the  meteorological  condi¬ 
tions  covered  so  far  in  flight  tests  are  limited  arid  do  not  include  measurements 
under  strong  convective  clouds  or  at  very  low  altitudes  where  conditions  may  give 
rise  to  a  somewhat  different  spectral  form. 

3.2.3  Stationarity  and  Homogeneity 

As  indicated  earlier,  properties  of  stationarity  and  homogeneity  of  a  random 
process  specify  an  invariance  in  the  statistical  characteristics  of  turbulence  with 
both  starting  time  and  starting  position,  respectively.  Inasmuch  as  the  intensity 
of  turbulence  and,  in  fact,  its  very  existence  in  the  atmosphere  is  dependent  upon 
the  broader  weather  conditions,  these  properties  can  only  apply  in  a  very  limited 
sense.  Weather  conditions  generally  involve  large-scale  organized  patterns  of  air 
motions  which  frequently  cover  hundreds  of  miles  and  which  change  relatively  slowly 
with  time;  significant  changes  ordinarily  take  place  in  a  matter  of  hours.  As  a 
consequence,  homogeneous  and  stationarity  conditions  might  at  best  only  be  expected 
to  apply  within  limited  regions  of  perhaps  the  order  of  100  miles  and  tii?*  durations 
of  perhaps  the  order  of  1  hour.  Some  of  the  available  data  that  have  a  bearing  on 
tne;  questions  are  presented  as  an  indication  of  the  applicability  of  tu?*e  properties. 

• 

In  Figure  6,  a  comparison  is  given  of  the  power  spectra  of  vertical  gust  velocity 
obtained  from  two  successive  sections  of  a  flight  record.  These  individual  spectra 
each  covered  a  flight  distance  of  about  15  miles.  The  total  record  time  covered  by 
the  two  records  is  roughly  4  minutes.  The  frequency  range  covered  in  these  tests 
is  quite  limited  but  covers  the  gust  wavelengths  of  about  300  to  6,000  ft,  which 
are  the  more  important  gust  wavelengths  for  most  airplane  turbulence  response 
problems.  The  two  spectra  agree  in  general  form,  considering  that  there  are 
differences  in  both  time  and  space  involved  in  these  records,  However,  ths  spectra 
also  display  a  significant  variation  in  the  intensity  of  the  turbulence  at  the  lower 
frequencies,  with  the  turbulence  of  pert  I,  overall,  being  mere  severe  than  that 
for  part  II.  One  might  surmise  from  this  result  that  there  is  at  least  some  degree 
of  both  homogeneity  and  stationarity  present  in  the  atmosphere  at  the  higher 
frequencies,  although  variations  do  appear  to  exist  even  locally  for  the  lower 
frequencies. 

3.2.4  Gauss  ianness 

The  next  property  to  be  considered  is  that  of  the  Gaussianness  of  atmospheric 
turbulence.  As  indicated  in  the  Appendix,  the  condition  of  Gaussianness  implies 
that  the  disturbance  fluctuations  have  a  Gaussian  probability  distribution. 

Furthermore,  a  Gaussian  random  process  is  one  in  which  the  joint  distribution  of 
the  velocity  fluctuations  at  several  (all)  points  has  a  joint  Gaussian  distribution. 

The  rigorous  testing  of  experimental  data  for  these  conditions  is  by  no  means  a 
simple  undertaking.  For  present  purposes,  it  will  suffice  to  consider  merely  the 
overall  probability  distributions  of  the  velocity  fluctuations.  This  limited 
consideration  is  by  no  means  conclusive  but  is  generally  a  satisfactory  practical 
guide. 
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A  number  of  such  evaluations  have  been  made  in  the  literature  (see,  for  example. 
Reference  23).  Some  recent  results  obtained  by  the  N.A.C.A.  on  this  question  are 
representative  of  these  results  and  are  shown  in  Figure  7.  The  figure  shows  the 
cumulative  frequency  distribution  of  gust  velocity  plotted  on  probability  paper. 

The  paper  is  scaled  in  a  manner  that  yields  a  straight  line  for  a  normal  or 
Gaussian  probability  distribution.  The  data  plotted  were  based  on  1/10-sec  readings 
for  a  4-min  test  run.  For  comparison  with  the  measured  distribution,  a  fitted 
normal  distribution  is  also  shown.  It  is  readily  seen  that  the  fitted  curve 
approximates  the  measured  distribution  quite  well  although  some  discrepancy  is 
evident.  Other  comparisons  of  the  same  type  indicate  similar  results  although,  on 
some  occasions,  the  measured  distributions  appear  to  depart  significantly  for  the 
extreme  values  from  what  would  be  expected  for  a  Gaussian  distribution.  However, 
for  many  practical  purposes,  the  Gaussian  approximation  appears  quite  adequate  and 
reasonable. 

3.2.5  Isotropy 

As  mentioned  earlier,  isotropy  places  a  number  of  important  restrictions  on  the 
allowable  turbulence  variations.  The  first  property  associated  with  isotropy  is  the 
property  of  the  invariance  of  the  statistical  characteristics  of  turbulence  with 
airplane  flight  direction.  Figure  8  shows  some  results  obtained  in  Reference  11  on 
the  spectrum  of  vertical  gust  velocity  as  measured  in  two  successive  flight  runs  in 
an  upwind  and  a  crosswind  flight  path.  These  measurements  unfortunately  do  not 
cover  the  whole  frequency  region  of  interest  but  are  restricted  to  the  higher 
frequencies  covering  gust  wave  lengths  of  10  ft  to  several  hundred  ft.  However,  for 
this  region  of  gust  wavelengths,  the  spectra  are  in  very  good  agreement,  suggesting 
that,  at  least  for  the  frequencies  covered,  the  condition  of  isotropy  was  well 
approximated. 

A  second  property  of  isotropy  turbulence  which  is  important  is  the  specific 
relations  required  between  the  statistical  characteristics  of  various  turbulence 
components.  As  an  illustration  of  the  characteristics  required,  Figure  9  shows 
the  lateral  (or  vertical)  power  spectrum  of  gust  velocity  defined  by  Equation  10 
as  well  as  the  associated  spectrum  for  the  longitudinal  gust  velocity  specified  by 
Equation  (24)  in  the  Appendix  for  the  use  of  isotropic  turbulence.  From  Appendix 
Equation  (24),  it  may  be  shown  that  the  longitudinal  component  designated  in  the 
Appendix  oy  F(ft)  is,  for  this  case,  given  by 

; — w? 

v  v  v  1  +  QzLz 

For  isotropic  turbulence,  both  the  vertical  and  Side  component  of  the  turbulence 
sensed  by  an  airplane  would  be  expected  to  have  the  same  spectrum.  The  longitudinal 
or  head-on  component  of  turbulence,  while  having  the  same  mean-square-gusfc  velocity, 
does,  however,  have  a  different  spectral  form  such  as  indicated  for  the  sjeeial  case  by 
the  foregoing  expression  and  as  illustrated  in  Figure  9.  As  a  consequence  of 
Equation  (24)  of  the  Appendix,  measurements  of  a  single  component  of  turbulence  can 
be  used  to  determine  the  spectra  of  the  other  components  for  the  case  of  isotropic 
turbulence.  In  addition,  for  the  case  of  isotropic  turbulence,  the  measurement  of  a 
single  power  spectrum  also  permits  the  determination  of  cross  spectra  between 
directional  components  at  any  two  points.  The  required  relations  may  be  obtained 
from  Equation  (22)  of  the  Appendix  and  are,  for  example,  given  in  Reference  26. 
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A  few  experimental  results  on  the  question  of  isotropy  are  shown  in  Figures  10 
and  11.  In  Figure  10.  a  comparison  is  shown  of  the  power  spectra  of  the  vertical 
and  the  side  (or  lateral)  components  of  turbulence  as  obtained  in  Reference  11. 

The  agreement  between  these  two  power  spectra  is  quite  good  over  the  frequency 
region  covered  although  the  lateral  or  side  component  appears  more  severe  than  the 
vertical  component  of  the  turbulence.  As  in  the  previous  case,  the  full  region  of 
frequencies  of  interest  is  not  covered,  with  the  data  being  restricted  to  the  higher 
frequencies  only.  Some  additional  data  recently  obtained  at  the  N.A.C.A.  are  shown 
in  Figure  11.  The  curves  shown  here  represent  the  spectra  obtained  for  the  vertical 
and  the  longitudinal  turbulence  components.  The  longitudinal  spectrum  was  obtained 
from  measurements  of  a  sensitive  and  high-frequency  airspeed  system.  Corrections 
for  airplane  longitudinal  acceleration  were  not  made. ,  As  indicated  earlier,  for 
isotropic  turbulence,  che  spectrum  of  the  vertical  and  longitudinal  components 
should  differ  in  the  manner  indicated  in  Figure  9.  At  the  higher  frequencies,  there 
appears  little  difference  between  the  two  spectra  of  Figure  11.  At  the  'lower 
frequencies,  the  differences  appear  large,  but  these  results  are  probably  not  very 
reliable  in  this  frequency  region. 

As  a  matter  of  interest,  it  might  be  mentioned  that  efforts  were  also  made  to 
obtain  the  power  spectrum  of  the  lateral  (or  side)  gust  component  for  the  test  data 
covered  in  Figure  11.  Preliminary  results  obtained  appeared  to  indicate  that  at  the 
lower  frequencies  the  lateral  turbulence  was  substantially  more  severe  than  the 
vertical.  The  reliability  of  this  set  of  measurements  has,  however,  not  yet  been 
established. 


It  appears  from  the  foregoing  that  at  best,  the  condition  of  isotropy  nay  be 
expected  to  apply  to  atmospheric  turbulence  to  a  limited  degree.  In  particular, 
the  available  Information  suggests  that  isotropy  nay  only  be  approximated  at  the 
higher  frequencies.  At  the  lower  frequencies,  anisotropic  conditions  apnear  to  be 
indicated  by  at  least  some  of  the  measurements  available. 

3.2.6  Limitations 

Mention  should  be  made  of  the  limitations  present  in  the  data  so  far  obtained. 

Fox-  the  most  part,  the  turbulence  measurements  have  been  restricted  to  limited 
frequency  regions  and  to  a  few  atmospheric  conditions.  Measurements  were,  in  most 
cases,  obtained  under  clear-air  flight  conditions  at  moderate  altitudes  from 
1,000  to  5,000  ft  above  terrain.  There  is  good  reason  to  believe  that  under  other 
atmospheric  conditions,  the  turbulence  characteristics  might  be  quite  different  from 
those  so  far  obtained.  First,  for  very  low  altitude  conditions,  say  below  1,000  ft, 
it  might  well  be  expected  that  turbulence,  being  frequently  wind  generated,  would 
more  closely  reflect  the  local  terrain  characteristics.  As  a  consequence,  conditions 
of  homogeneity  would  be  largely  dependent  upon  terrain  homogeneity  and  conditions  of 
stationarity  would  be  dependent  upon  persistence  of  the  broad  wind  field.  The 
spectral  characteristics  for  this  case  might  also  be  expected  to  be  more  sensitive 
to  terrain  as  well  as  the  thermodynamic  conditions.  Finally,  it  is  questionable 
whether  conditions  of  isotropy  might  be  approximated  under  such  conditions  and.  in 
fact,  some  measurements  now  available  fre®  meteorological  towers  suggest  that 
isotropy,  while  approxinated  at  the  higher  frequencies,  does  not  appear  to  extend 

over  the  entire  region  of  frequencies  that  is  of  interest  in  aeronautical  applica¬ 
tions. 


A 
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The  second  weather  condition  that  deserves  special  comment  is  that  of  strong 
convective  clouds  or  thunderstorms.  Because  of  the  severe  conditions  of  turbulence 
under  such  conditions,  this  weather  condition  is  of  particular  concern  in  aeronautical 
applications.  The  basic  physical  element  of  a  thunderstorm  is  normally  of  small 
dimens^or.,  of  the  order  of  perhaps  5  to  10  miles  in  height  and  width,  and  contains 
a  we II -organized  large-scale  flow  pattern  which  generates  smaller  eddies  or  turbulent 
velocity  fluctuations29.  As  a  consequence,  it  might  be  expected  that  the  turbulence 
would  more  closely  reflect  the  larger  scale  thunderstorm  structure  and  turbulence 
measurements  obtained  under  ether  conditions,  such  as  in  clear  air,  may  not  be 
representative  of  conditions  for  this  case. 


4.  PROBABILITY  DISTRIBUTIONS  OF  ROOT -MEAN -SQUARE  GUST  VELOCITY 

In  the  preceding  section,  it  was  indicated  that  for  given  conditions  of  time  and 
place,  atmospheric  turbulence  might  be  considered  to  approximate  a  homogeneous  and 
Gaussian  random  process.  Over  mest  of  the  frequency  region  of  concern,  the  spectral 
form  of  the  vertical  (or  lateral)  gust  velocity  appeared  to  be  approximated  by  the 
following  expression  (Equ.  10) 


$<fl)  = 


_  2  _L  1  +  3  02L2 
w  v  (i  +  d2l2)2 


with  values  of  the  scale  of  turbulence  L  of  the  order  of  1,000  ft.  However,  the 
values  of  the  root-mean-square  gust  velocity  cr#  appear  to  vary  widely  with  weather 
condition.  In  this  section  of  the  present  paper,  this  simplified  model  of  the  local 
turbulence  will  be  used  as  a  building  block  to  provide  a  more  general  description 
of  the  overall  turbulence  history  experienced  by  airplanes  in  operational  flight. 
Such  a  description  is  desirable  and  needed  in  order  to  permit  the  determination  of 
the  airplane  loads  and  other  response  histories  in  operational  flight. 


The  detailed  measurements  of  turbulence  and  turbulence  spectra  discussed  in  the 
last  section  are  still  far  too  meager  to  supply  the  information  required.  Direct 
measurements  of  turbulence  power  spectra  to  provide  what  might  be  called  a  turbulence 
climatology  appears  to  be  a  distant  future  possibility.  A  vast  body  of  statistical 
data  from  normal  operations  on  airplane  acceleration  responses  in  rough  air  does 
exist,  however.  These  data  have,  in  most  cases,  been  collected  in  terms  of  the 
total  number  of  peak  accelerations  (or  effective  gust  velocities)  that  exceeded 
given  values  of  acceleration  (or  gust  velocity)  in  operational  flights.  Under 
certain  simplifying  assumptions,  it  has  appeared  possible  to  use  these  data  to 
derive  reasonable  estimates  of  the  turbulence  climatology  of  the  atmosphere  in  the 
form  of  the  probability  distributions  of  the  root -mean-square  vertical  gust  velocity 
for  various  conditions.  Such  results  are  reported  in  References  18,  19  and  25. 
Inasmuch  as  these  results  provide  at  leant  a  beginning  of  a  description  of  tirbulence 
that  can  be  used  in  spectral-type  calculations  of  operational  response  histories,  it 
appeared  appropriate  to  review  briefly  the  general  approach  used  and  some  of  the 
results  obtained. 


I 
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Method  of  Analysis 

» 

The  general  nethod  used  In  theee  studies  Involves  t»o  steps.  These  are: 

(1)  the  derivation  of  the  probability  distribution  of  the  root-mean-square 
accelerations  from  the  data  on  counts  of  acceleration  peaks 

(2)  the  transformation  of  the  root -mean -square  acceleration  to  root -mean-square 
gust  velocities  through  the  use  of  the  airplane  response  characteristics. 

These  two  steps  will  be  discussed  in  order. 

k.l.l  Derivation  of  the  Probability  Distribution  of  Root -Mean-Square 
Acce  leration 

For  a  stationary  Gaussian  random  process,  say  y(t),  a  simple  approxinate  relation 
exists  between  the  number  of  peaks  per  second  exceeding  a  given  value  of  y  and  the 

power  spectrum  of  y  and  is  derived  in  Reference  26.  This  relation  is  also  discussed 
in  References  18  and  25  and  is  given  by  discussed 

N(y)  -  N'  e-y2/2aJ 

-  "oe  (11) 

where  a2  equal  to  the  mean  square  value  of  y  and  N  can  be  considered  a  character¬ 
istic  frequency  of  y(t)  and  is  defined  by  cnaracter 


N0  = 


1 

“  00  " 

I  oo2$  (co)dco 
0  y 

2  TT 

QD 

/  ^v(oi)doj 

-  o  J 

(11a) 


In  jteference  18,  this  relation  is  extended  to  the  case  of  an  airplane  flying 

W  ™i™w””SrrS£  ,lr  °f  "*ed  SMCtrl1  f0r"  but  r°ot-mean-square 

exceeding  “given ““''"“h10"8  *»  - 


M(a„)  =  N0  f 


where 


P(cr_)e 


2~ 


(12) 


M(an>  =  ^IduTo/T  nUmb6r  °f  1)64)1  acceleratlons  Per  second  exceeding  a  given 

°a  =  root -mean-square  value  of  an 

P(CV  =  the  probability  distribution  of  cr  . 

a 

Equation  12  permits  the  determination  of  the  probability  distribution  of  th* 
mean-square  acceleration  p(cr  )  from  measured  distribution  of  acceleration  peaks  1(C) 
in  most  cases,  the  Integral  equation  (13)  cannot  be  solved  In  closed  fore  ^ut  % 

ms  .I"  6StI,KteS  °f  P<<V  *  "-»1»  derived °(see ^Reference  25) 

(This  same  approach  may  also  be  applied  if  the  Measurements  are  In 

S:S‘iM  *  °0nSl,ierln8  tte  'ftbbt!"  -'-‘ties  as  reduced 


4 .1.2  Probability  Distributions  of  Root -Mean-Square  Gust  Velocity 

The  conversion  of  the  derived  distributions  of  root -mean-square  acceleration  to 
those  for  root -mean-square  gust  velocity  requires  the  determination  of  the 
relation  between  the  two  root -mean-square  values.  If  the  airplane  is  assumed  a 
linear  system  and  average  values  of  such  airplane  characteristics  as  weight,  air¬ 
speed,  air-density,  center-of -gravity  position,  etc.,  are  taken  as  adequate,  then 
for  a  given  operation  or  portion  of  an  operation,  the  root -•mean-square  values  are 
related  simply  as  follows 


where 


(13) 

l 

?  (14) 


For  the  case  of  a  rigid  airplane  free  to  move  vertically  but  not  pitch,  the  quantity 
A  is  given  by 


A  -  ■PVStn  /i(K.s) 
2W  V  rr 


where 


p  -  air  density 

V  =  forward  speed 


S  =  wing  area 

m  =  slope  of  the  lift  curve 


* 


i) 


W  =  airplane  weight 

I(K,s)  =  an  airplane  gust  response  function  (defined  in  Reference  4)  which  depends 
upon  the  airplane  mass  parameters  =  4W/gnpSc  and  s  =  c/L  where  c  is  the 
mean  wing  chord. 


This  case  is  used  in  References  18,  19,  and  25. 

If  the  relation  between  root -mean-square  values  of  Equation  (13)  is  assumed  to 
apply,  the  probability  distribution  of  root -mean-square  gust  velocity  is  then  given 
by 

P(crw)  =  Ap(ir&)  (16) 

Equation  (16)  permits  the  simple  conversion  of  probability  distributions  of  root- 
mean-square  accelerations  to  root -mean-square  gust  velocities.  This  conversion 
step  is,  for  the  case  considered,  seen  to  depend  only  upon  the. response  quantity  A. 


4*2  Application  to  Guat  Load  Data 
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In  this  section,  some  results  obtained  by  the  aDDliratinn  ,  .... 

In  the  preceding  section  to  sv.ll.Dle  statistic!  "t, « 1  ,!£  r°1,tl0ns  derl,ed 
Three  principal  problem  are  encountered  ln\«“  a™  wif™  ™  T  ‘"'eSe"ted- 
•choice  of  an  appropriate  value  of  the  scale  « Z Z T»  lo7  Z7  “* 
min. t Ion  of  the  airplane  acceleration  freuuency  respo^  fnneH™  ?  '  .  Stei" 

more  particularly,  the  value  of  A)  and  tha  fl^tlon  for  gusts  H(oj)  (or 

,  _  '•  the  determination  of  the  val'je  of  th»  ehomo 

terlstic  frequency  N  .  The  value  nf  r  _  *  *»Aje  or  me  cnsrac- 

Cm  ithmiffh  ,  0  ue  of  L  chosen  for  present  purposes  was  1  000  ft 

(although  other  values  were  also  considered  in  Reference  181  whioh 

'  r  1~“?  r 

;“fht°  ^ ~ 

for  the  various  airplane.  In  .chance  siVtT 


4  .2.1  Variation  of  p(aw)  with  Altitude 


In  order  to  arrive  at  some  estimates  of  the  variation  of  6 (a  \  with  „ 

was  made  of  the  sunm»ry  of  aist  .  °  0  P(  w>  lth  altltude.  use 

Be  fevers  „  ,  8U™*ry  or  gust  statistics  given  in  Reference  31.  Figure  6  of 

shich  are  In  term,  of  derived  gust  velocities  .ere  hT™.?  o  tS' 

£*££"*£,*  ““  “^‘"istlcs'of  .  representative  "t  raw  port  Z'laT 

^-rss?iLrrMit*ar  ' 

altitude  bracket,  the  number  of  peak  accelerations!  for  the  various  altltnto  >7.  , 
could  be  approtlmnted  by  the  follwlng  distributions  for  6<o>):  blanket. 

for  the  altitude  bracket  of  0  to  10,000  ft: 


P(a,)  =  0.99-1.  e^w/l,Ha+  Q.oi  _J_  e'0"72-8* 


1.48 


2.48 


(17) 


for  the  altitude  bracket  10,000  to  30.000  ft: 


£(ct  )  =  - 1 -  e"1^*^0,32 

w  2(0.32) 2 

for  the  altitude  bracket  40,000  to  50,000  ft: 


(18) 


P(e.)  = 


~V&I/0.  2  9 


2(0.29)2 


(19) 
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These  distributions  are  shown  in  Figure  12  in  terms  of  their  cumulative  probabil 
ity  distributions  P (crw)  (obtained  by  integrating  Eqs.  17  to  19  from  given  values  of 
crw  to  infinity)  and  indicate  the  proportion  of  flight  time  spent  above  given  valties 
of  o'w  for  the  various  altitude  brackets.  Perhaps  the  most  important  points  to  be 
noted  in  Figure  12  are  the  relatively  large  amount  of  time  spent  in  essentially 
smooth  air  at  the  higher  altitudes  (c^  <  2  ft/sec,  93  to  96%  of  the  time)  and  the 
relatively  large  amount  of  time  spent  in  light  to  severe  turbulence  at  the  lowest 
altitude  bracket  (crw  >  2  ft/sec,  25%  of  the  time).  The  total  time  spent  for 
ctw  >  5  ft/sec  for  the  lowest  altitude  bracket  is  roughly  5  to  iC  times  as  great 
as  that  for  the  higher  altitude  brackets. 


4.2.2  Method  of  Application 

o 

The  foregoing  results  may  be  applied  in  a  straightforward  manner  to  the  calcula¬ 
tion  of  response  histories  for  arbitrary  operations.  The  method  of  application  of 
these  results  to  the  calculations  of  response  histories  is  based  on  the  following 
relation: 


dcrw  (20) 

which  is  obtained  by  substituting  Equations  13  and  16  into  Equation  12.  The  procedure 
involves  the  division  of  the  operational  flight  plan  into  homogeneous  portions  or 
segments  in  regard  to  flight  altitude  and  operating  conditions  such  as  airspeed  and 
airplane  weight.  The  appropriate  distribution  of  the  root -mean-square  gust  velocity 
is  selected  for  each  flight  segment  from  Figure  12  and  Equation  20  is  then  evaluated 
for  each  segment.  Actual  numerical  calculations  are  facilitated  by  the  use  of  the 
charts  given  in  References  18  and  25.  The  sequence  of  steps  involved  in  soch 
applications  is  as  folltwrs: 


“<v  -- 


A  -a.  72  (Ar  ) 

0(<Ve 


(1)  The  operational  flight  plan  is  divided  into  homogeneous  segments  in  regard  to* 
flight  altituje  (10.000  ft  altitude  bracket?)  and  operating  conditions  such  as 
airspeed  and  weight; 

a 

ft 

(2)  The  appropriate  distribution  of  p(aw)  is  selected  for  each  altitude  bracket 
from  Figure  12; 

(3)  The  values  of  A  for  the  acceleration  or  other  response  are  determined  for 
each  significant  segment  of  the  flight  plan  in  accordance  with  Equation  (14); 

(4)  In  order  to  obtain  the  associated  distributions  of  loot -mean-square  accelera¬ 
tion  P(°a),  each  of  the  distributions  of  p(crw)  is  transfor>«ed  by  the  relation 

p(cra)  =  yP(o"w) 
cr 

where  crw  =  _JL 

A 

(5)  The  values  of  NQ  are  most  easily  determined  from  flight  records,  if  available, 
by  the  methods  already  indicated25.  For  new  designs  N0  must  be  estimated 
analytically  by  application  of  Equation  (11a).  In  terms  of  the  power  spectrum 
of  the  gust  input,  NQ  is  given  by 
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(6)  The  distribution  of  and  the  values  of  NQ  are  then  used  in  Equation  (12)  to 
derive  the  number  of  peak  accelerations  per  second  or  per  mile  for  each 
condition.  These  calculations  are  facilitated  by  the  use  of  charts  such  as 
given  in  Reference  (18)  or  (25). 

(7)  The  results  obtained  in  step  6  are  then  weighted  in -accordance  with  the  flight 
distance  in  each  condition  or  segment  and  then  summed  for  all  conditions  in 
order  to  obtain  the  overall  acceleration  or  other  response  history. 


5.  CONCLUDING  REMARKS 


The  foregoing  discussion  has  served  to  provide  a  preliminary  assessment  of  the 
characteristics  of  continuous  atmospheric  turbulence  as  determined  from  recent 
flight-test  measurements.  In  the  discussion,  the  general  methods  that  have  been 
used  for  obtaining  measurements  of  continuous  atmospheric  turbulence  have  been 
reviewed  and  some  of  their  limitations  indicated.  In  spite  of  the  limitations,  the 
techniques  used,  and  particularly  the  technique  based  on  flow  direction  measurements, 
appear  to  provide  a  reliable  measurement  of  continuous  atmospheric  turbulence.  In 
evaluating  the  measurements  so  far  obtained,  the  local  characteristics  of  turbulence 
were  considered,  and  particular  attention  was  given  to  the  properties  of  the  power 
spectrum  of  turbulence.  PGr  the  measurements  so  far  obtained,  the  power  spectrum 
appears  to  show  a  persistent  spectral  form  over,  a  wide  range  of  frequencies  and,  in 
particular,  appears  to  decrease  with  increasing  frequency  as  l/Q 2  for  the  range  of 
frequencies  of  principal  concern. 


The  following  analytic  expression 


rr  2_L  1  +3H2L2 
W  n  (1  +  D2L2)2 


where 

ctw  =  root -mean-square  gust  velocity,  and 
L  =  scale  of  turbulence, 

appears  to  represent  the  available  measurements  of  the  spectrum  of  vertical  gust- 
velocity  reasonably  well.  For  the  available  measurements,  the  value  of  L  appears 
to  be  about  1,000  ft,  while  values  of  a  vary  over  a  wide  range. 

*  W  * 

In  addition  to  the  spectral  properties,  the  characteristics  of  the  measurements 
were  also  examined  in  regard  to  properties  frequently  used  in  idealised  and  simplified 
theoretical  models.  In  particular,  the  properties  of  stationarity,  homogeneity, 
Gaussianness,  and  isotropy  were  considered  and  the  measurements  examined  to  determine 
how  well  these  properties  apply  to  the  atmosphere.  Examination  of  the  measurements 
available  suggests  that  these  properties  appear  to  apply  to  atmospheric  turbulence 
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to  a  limited  degree.  The  applicability  of  most  of  these  properties  appears  best  at 
the  highest  frequencies  and  poorest  at  the  lower  frequencies.  The  evidence  thus 
suggests  that  while  analytical  studies  based  on  such  simple  uodels  may  prove  useful 
guides,  they  must  be  applied  with  caution.  Further  Measurements  of  atmospheric 
turbulence  are  needed  in  order  to  cover  both  a  wider  range  of  frequencies  and  a 
wider  range  of  meteorological  conditions.  In  particular,  a  need  for  measurements 
under  such  important  conditions  as  very  low  altitudes  and  strong  convective  clouds 
or  thunderstorms  is  indicated. 

In  addition  to  the  consideration  of  the  local  characteristics  of  atmospheric 
turbulence,  consideration  was  also  given  to  the  problem  of  applying  the  techniques 
of  continuous  turbulence  analysis  to  the  calculation  of  airplane  response  in 
operational  flight.  For  this  purpose,  recent  results  obtained  on  the  probability 
distributions  of  the  root -mean-square  gust  velocities  were  described  and  their 
method  of  application  to  response  calculations  indicated. 


♦ 
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Pig.  1  Method  of  detemination  of  the  power  spectrum  of  gust  velocity  from  airplane 

response  measurements 


(3)  wg=  V(av-fl)-w0  +  /d 


« 

Pig.  2  Geometric  relations  between  vertical  gust  velocity  and  vane  *ngle 


ft 


Pig.  3 


CMparlsm,  of  peer  spectr*  of  gust  velocity  as  obtained  ,1th  and  without 
corrections  for  airplane  pitching  and  vertical  velocity  notions 


«»,  FT/ SEC 


V,  FT/SEC 


0n,  FT/SEC2 


Pig. 4  Comparison  of  time  history  of  vertical  gust  velocity,  side  gust  velocity, 

and  airplane  nornail  acceleration- 
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Pig.  5 


Fig. 
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Measured  power  spectrum  of  vertical  gust  velocity  and  analytic  representation 


fl,  RAO/FT 

Comparison  of  power  spectra  of  vertical  gust  velocity  for  two  parts  of  a  test  run 
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« 


Pig.  8 


Comparison  of  the  power  spectra  of  vertical  gust  velocity  for  upwind  and 

crosswind  flight  surveys  (Ref.  11) 
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Pig. 9  Comparison  of  lateral  (or  vertical  gust  velocity  power  spectrum  of  Equation  10 
in  the  main  text  with  longitudinal  spectrum  for  case  of  isotropic  turbulence 


Pig.  10  Comparison  of  the  power  spectra  of  lateral  and  vertical  gust  velocity  (Ref.  11) 

*  *  * 
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Pig.  11  Comparison  of  the  power  spectra  of  vertical  and  longitudinal  gust  velocity 


Fig. 12  Cumulative  probability  distributions  of  'oat- mean-square  gust  velocity  for 

various  sltitude^ 
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APPENDIX 


Fundamental  Concepts  of  Random  Process  Theory 


1.  GENERAL  REMARKS 

In  this  appendix,  some  of  the  basic  concepts  and  relations  of  random  process  theory 
and  the  statistical  theories  of  turbulence  pertinent  to  the  analysis  and  interpretation 
of  measurements  of  atmospheric  turbulence  will  be  reviewed.  Rather  complete  treatments 
of  these  two  areas,  random  process  theory  and  the  statistical  theories  of  turbulence, 
are  contained  in  References  30  and  27,  respectively.  The  present  discussion  is  aimed 
at  defining  some  of  the  specialized  terminology  and  indicating  a  few  of  the  more 
important  relations  that  are  required  for  present  purposes.  The  discussion  will  start 
with  the  consideration  of  a  simple  random  process  f(x)  involving  the  single  parameter 
x.  The  basic  concept  of  a  power  spectrum  of  a  random  process  and  the  important  case 
of  the  Gaussian  random  process  will  be  described.  An  indication  will  also  be  given 
of  the  methods  of  application  of  these  ideas  to  airplane  response  problems. 

In  many  aeronautical  applications,  generalizations  of  the  foregoing  simple  cases 
are  required.  For  example,  recent  studies12'22  have  indicated  that  in  addition  to 
variations  oil  turbulence  along  the  flight  path,  the  variations  in  turbulence  in  the 
airplane  spaawise  direction  must  frequently  be  taken  into  account  in  determining 
both  airplane  lateral  and  longitudinal  responses.  The  consideration  of  the  turbulence 
variations  in  these  two  directions  requires  an  extension  to  two-dimensional  processes 
f(x,y).  The  most  general  case  that  might  bo  considered  is  a  four-dimensional  process 
f(x,y,z, t);  the  fourrdimensional  process  accounting  for  the  variations  in  all  three 
space  coordinates  as  well  as  time.  However,  in  practice,  this  full  generalization 
is  rarely  required  and  would  unduly  complicate  the  present  discussion.  Turbulence 
variations  with  time  can  frequently  be  neglected  on  the  basis  that  the  turbulence 
pattern  changes  slowly  with  respect  to  the  rapid  forward  speed  of  the  airplane. 

This  assumption  is  the  equivalent  of  Taylor's  hypothesis,  which  is  commonly  applied 
in  studies  of  wind-tunnel  turbulence.  Thus,  in  aeronautical  application,  the  air¬ 
plane  is  generally  envisioned  as  sampling  an  instantaneous  cross  section  of  the 
turbulence  in  space;  The  variations  of  turbulence  in  the  z  direction  mey  also  usually 
be  ignored  since  in  most  cases  the  airplane  motions  in  a  direction  normal  to  the 
average  flight  path  are  considered  small  * 

In  addition  to  the  generalization  to  multidimensional  random  processes,  a  further 
generalization  to  vector  processes  is  frequently  required  in  order  to  account  for  the 
directional  components  of  tho  turbulence.  For  example,  in  regard  to  the  airplane 
normal  response,  an  adequate  analysis  sometimes  requires  the  consideration  of  both  the 
normal  and  the  longitudinal  or  head-on  components  of  atmospheric  turbulence.  Thus, 
in  the  general  case,  the  process  f(x,y)  has  to  be  considered  as  a  vector  process  with 
three  normal  directional  components  u,  v,  and  w.  The  general  consideration  of  multi¬ 
dimensional  vector  processes  rapidly  becomes  too  involved  for  many  practical  appli¬ 
cations.  As  a  consequence,  the  necessity  exists  for  simplifications  if  practical 
results  are  to  be  achieved.  Fortunately,  simplifications  can  frequently  be  made 
involving  symmetries  in  the  space  coordinates  and  internal  relationships  between  the 
statistical  characteristics  of  the  turbulence  directional  components.  The  most 
important  and  useful  of  these  simplifications  is  the  case  of  isotropic  turbulence, 
which  will  also  be  considered. 


2.  SINGLE  DIMENSIONAL  RANDOM  PROCESSES 


2.1  Random  Processes 

The  term  random  immediately  implies  a  lack  of  definiteness  that  can  only  be  des¬ 
cribed  in  statistical  or  probability  terms.  In  the  general  sense,  a  random  process  * 
may  be  said  to  consist  of  a  set  of  functions  (ensemble),,  say  f(x),  where  the  para¬ 
meter  x  may,  for  example,  designate  time  or  as  in  the  present  case  a  space  coordinate, 
and  where  the  individual  functions  fa(x),  fb(x),  etc.,  each  occurs  with  a  given 
probability.  Alternately,  the  process  f(x)  may  be  defined  in  terms  of  probability 
distributions  for  given  values  of  x.  The  function  p(fx)  is  used  to  designate  the 
probability  distribution  of  f(x)  at  position  x.  and  the  joint  probability  distribu¬ 
tion  of  f  for  various  positions,  x1#  x2,  etc. ,  ^v  be  designated  by 

P(fx,  f2,  .  .  .) 

where  the  subscripts  designate  the  values  of  f  for  various  values  of  the  parameter  x. 
The  probability  distribution  may  be  said  to  define  the  proportion  of  all  members  of 
this  set  that  have  given  combinations  of  values  at  the  various  positions,  and  is  thus 
taken  to  define  the  statistical  characteristics  of  the  process  completely. 


The  definition  ot  a  random  process  in  terms  of  its  complete  probability  distribu¬ 
tion  is,  however,  generally  too  cumbersome  for  practical  applications  and  recourse 
is  generally  had  to  simpler  and  more  manageable  properties  of  the  process.  Of  these, 
the  most  important  are  the  average  values 

Av.  {f (x) }  =  /*  fp(f)df 

-  00 


Av.  {f2(x)>  =  /*  f2p(f)df 
-  00 

Av.{f(xl)  f(x2)>  =  /eBf1f2p(f1,t2)df1  df2 


V 

y  (a, i) 


=  Rff<*l»*2>  J 

These  three  averages  are  sufficient  and  most  useful  for  niauy  purposes.  The  first 
expression  defines  the  average  value  for  various  values  of  x.  The  second  quantity, 
Av.  {f2(x)},  defines  the  mean-square  value,  and  the  quantity  Av. {f(Xj)  f(x2)},  defines 
the  average  product  of  ft  and  f2.  In  the  special  case  of  a  stationary  Gaussian 
random  process,  which  is  of  particular  interest  in  present  considerations,  these 
quantities  will  be  seen  to  specify  the  statistical  characteristics  of  the  process 
completely. 


2.2  Stationarity  and  Homogeneity 

As  defined  in  Equation  (A. 1)  herewith,  the  average  values  are  clearly  functions 
of  the  parameter  x.  In  many  practical  applications,  it  is  useful  to  consider  the 
special  case  in  which  these  averages  are  independent  of  the  value  of  x.  If  the 
quantity  x  designates  time,  such  processes  are  termed  stationary  end  the  average 
values  and  probability  distributions  do  not  depend  upon  the  starting  point  Xj. 


A-iv 


Alternatively,  if  the  parameter  is  a  space  parameter,  as  is  generally  the  case  in 
turbulence  applications,  the  process  is  termed  homogeneous.  Under  such  conditions 


Av.  {f(xi)}  -  Av.  {f  (Xj) } 

Av.{f2(xi)}  =  Av.  {f2(Xj)> 

Av.  (ffxj)  f(Xj)>  =  Rff(Xj  -  Xi) 

Under  certain  conditions  (which  generally  apply  in  practice),  stationary  processes 
have  the  important  and  useful  ergodic  property  of  the  equivalence  of  averages  over 
the  set  of  functions  with  averages  in  time  (or  in  space)  for  a  single  function. 
Thus  we  may  write  that 


Av. {f(x)}  =  lim 
X~*« 

Av. {fz(x)}  =  lim 

x-*® 


Av.  {f (Xj )  f(Xj)>  = 


t<x)ix 

I//™* 

lim  f(x)  f(x+Ax)dx 

x_®^y-x 


>  (A. 3) 

Rff(Ax) 


where  Ax  -  x.  -  Xj.  Thus,  a  single  member  of  the  set,  which  is  frequently  all  that 
.  ls  available  in  practice,  may  be  used  to  define  the  statistical  properties  of  the 
entire  seta  It  will,  without  loss  of  generality,  be  assumed  that  the  mean  value 
defined  by  Av.  {f(x)}  is  equal  to  zero.  Por  this  case,  the  second  moment  average 
Av.  \f  (x)}  is  generally  termed  the  variance  or  mean  square  or  power  and  the  quantity 
Rff(Ax)  is  termed  the  auto-covariance  function  or  the  auto-correlation  function. 

The  auto-correlation  function  is  symmetrical  about  the  origin  and  has  its  maximum 
value  at  Ax  =  0. 


2.3  Power  Spectra 

* 

The  auto-correlation  function  is  particularly  important  since  it  reflects  the 
frequency  characteristics  of  the  disturbance,  and  its  Fourier  transform  gives  rise”" 
to  the  power  spectrum  of  the  random  process.  The  power  spectrum  of  f(x)  is  defined 
by  the  reciprocal  relations 


$ff<n>  =  [Rff(Ax)e'in  AxJ  dA: 

Rff(Ax)  =  A^*^(0)ein  Ax_|df) 


>  (A. 4) 


where 


H  =  reduced  frequency,  radians/ft,  and 
Ax  =  space  displacement,  ft. 
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The  power  spectrum  is  a  real  and  positive  symmetrical  function  because  of  the 
symmetry  of  Rff(Ax)  and  is  a  particularly  useful  function  since  it  describes  the 
distributions  of  the  energy  of  f2(x)  with  frequency.  Because  of  the  symmetry  of 

and  Rff(Ax),  Equation  (A.  4)  herewith  may  be  simplified  to  reciprocal  cosine 
transformations. 

2.4  Cross  Spectra 

For  the  case  of  two  random  processes,  f(x)  and  g(x),  additional  averages  of  the 
form 


Av.  (f{x^)  g(xp}  =  RfgUj  -  Xj) 

are  also  of  interest.  Such  averages  reflect  the  correlation  between  the  two 
processes  and  thus  Rfg(Xj  -  Xj)  is  called  the  cross-correlation  function.  For 
stationary  and  ergodic  processes 

1  /X  * 

Av.  (fUj)  g(Xj)}  =  lim  2 xf  g(x  +  Ax>dAx  (A.5) 

X—  00  u  X 

The  cross  spectrum  is  not  symmetrical  and  is  easily  shown  to  have  the  property  that 


Rfg(  Ax)  =  Rgf(-Ax) 


(A.  6) 


The  Fourier  transform  of  Rffe.(Ax)  may  also  be  defined  and  yields  the  so-called  cross 
spectrum  $>fg(fi)  between  f(x)  and  g(x).  The  reciprocal  relations  between  Rfg(Ax) 
and  are  given  by 


Vn> 


d  A  x 


=  c(H) 


Rfg(Ax) 


1  r 

2 


iq(fi) 

$fg(n)eiQA; 


d  n 


>  (A.  7) 


j* 


The  cross  spectrum,  in  contrast  to  the  power  spectrum,  is  complex,  having  both  real 
and  imaginary  parts.  The  real  part  c(Q)  is  generally  termed  the  cospectrum  and 
provides  a  measure  of  the  in-phase  power  of  the  two  disturbances.  The  imaginary 
part  q(fl)  is  generally  termed  the  quadrature  power  and  provides  a  measure  of  the  90° 
out-of-phase  power  between  the  two  disturbances.  From  Equation  (A. 6),  it  may  be  shown 
that 


Vn>  =  V<*n> 


(A.  8) 


2.5  Input-CXitput-Relations 

The  power  spectrum  and  cross  spectrum  have  the  important  property  of  permitting 
linear  operations  in  simple  form.  The  response  of  an  airplane,  flying  along  the  x 
axis  at  a  given  speed,  to  an  arbitrary  gust  disturbance  f(x),  may  be  expressed  as 

* 
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z(x) 


(A.  9) 


=  J*f (X  -  Xj)  hCXjJdXj 


where  h(x)  is  the  response  of  the  airplane  to  a  unit  impulse  gust  at  x  -  0  and  is 

generally  termed  the  unit  impulse  response  function.  The  linear  relation  of 

Equation  (A. 9)  leads  to  the  following  relations  among  the  power  spectrum  of  the 

response  4>  (0)  and  the  cross  spectrum  between  the  disturbance  and  the  response 

zz 


=  $ff(n)iK<n)|2 

*tz®> 

=  $ff(0)H(H) 

where 

H<n> 

=  j;*[h(x)e-in*]dx 

(A.  10) 
(A.  11 ) 


and  is  generally  termed  the  frequency  response  function  and  defines  the  system 
response  to  unit  sinusoidal  gusts  of  frequency  Cl.  Thus,  it  may  be  seen  that  the 
integral  relation  of  Equation  (A, 9)  yields  simple  product  relations  between  the 
spectra.  These  simple  product  relations  for  linear  systems  provide  a  major  reason 
for  the  usefulness  of  spectral  methods. 


2.6  Gaussian  Random  Process 

The  special  case  of  stationary  random  process  that  is  of  particular  interest 
because  of  its  natural  simplicities  and  because  it  frequently  approximates  the 
conditions  of  natural  phenomena,  is  the  case  of  a  Gaussian  process.  In  the  present 
paper,  it  is  of  interest  to  see  how  well  this  model  approximates  the  condition  of 
atmospheric  turbulence.  This  special  case  has  the  property  that  the  probability 
distributions  is  Gaussian.  Specifically,  for  any  value  of  x 


P(f) 


1  p-f  z/2cr2 


(A. 12) 


where  a2  =  Rff(0).  The  joint  probability  distribution  for  several  values  of  x  also 
has  a  multi-variate  Gaussian  distribution.  Thus,  for  example,  the  joint  distribution 
of  f  for  Xj  and  x2  is  given  by 


P(f,,f,) 


27 t(ct“  -  R122) 


exp  - 


"V  -  2R12f-lf2  +CT‘f 


2t  2 
2 


2  (a11 


R12) 


(A. 13) 


where  R12  =  Rff(x,  -  xx). 

This  distribution,  as  well  as  higher-order  distributions  for  xx,  x2.x3,  etc.,  depends 
only  upon  the  second-order  averages  of  the  random  process  given  by  the  auto¬ 
covariance  function  R^^(Ax),  and  in  this  sense,  the  process  is  said  to  be  completely 
specified  by  the  auto-covariance  function  or  its  Fourier  transform,  the  power  spectral 
density  function. 


The  characteristics  of  a  Gaussian  process  have  been  investigated  extensively  and  a 
number  of  useful  relations  have  been  derived  for  this  type  of  process  which  permit 
the  estimation  of  statistical  quantities  of  frequent  interest  from  the  basic  quantity, 
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the  power  spectrum.  Many  of  these  relations  are  given  in  Reference  28,  and  a  few  of 
the  more  important  ones  in  aeronautical  applications  are  given  in  Reference  20. 

For  a  stationary  Gaussian  process,  the  ergodic  property  will  generally  apply,  and 
the  probability  distributions  obtained  from  a  single  function  should  also  be  Gaussian. 
Thus,  a  single  time  history  sample  can  be  examined  to  determine  if  the  process  is 
Gaussian. 


3.  MULTIDIMENSIONAL  AND  VECTOR  PROCESSES 


3.1  Two-Dimensiona !  Power  Spectra 

In  analogy  with  Equation  (A.  1)  of  this  Appendix  for  the  single-dimensional  pro¬ 
cess,  we  may,  for  two-dimensional  process  f(x,y)  note  that 


Av. {f (x, y) }  =  /  ®fp(f)df 

-00 

Av.  {f?(x,y)}  =  /°°  f  ?p(f)df 


(A. 14) 


Av.{f(x1,y1)  f(x2.y2)>  =  /” f jf 2p(f f2)df j  df2  J 

For  the  case  of  homogeneous  and  ergodic  processes,  the  property  of  equivalence  of 
space  and  probability  averages  generally  applies  and 


Av.  {f(x,,y)>  =  iim 


-  f  fX 

J-V  J. 


Y  J-X 


f (x,y)dx  dy 


•  i  fY  rx 

Av.  {f 2 (x, y) }  =  lim  —  | 

X-®  4XY  J-Y  J-x 

Y-oo 


f 2(x,y)dx  dy 


i  rY  rx 

f.  {f  (Xj.yj)  f (x2,y?)}  =  lim  —  f(x1>yi)  f(x,,y  )dx  dyt 

X— ®  J-Y  J-X 


►(A. 15) 


=  Rff(Ax,  Ay) 


where  Ax  =  x?  -  Xl  and  -Ay  =  y2  -  y y  J 

Two-dimensional  power  spectra  may  also  be  defined  by  the  reciprocal  relations: 

$ff(fli,n2)  =  “5  [*  [  [nffiAx,  Ay)e'i^x+^2^y'  dAx  dAy 

J-GD  J»OoL  J 


Rff(  Ax,  Ay) 


» 

O0D 

-00 


n2)ei(niAx+nA)  dnv  dn2 


(A. 16) 
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3.2  Input-CXitput  Relations 


turbulence &f  (x  the ***  ^  X  &Xl<?  ^  &  giVe°  Spe8d  through  two-dimensional 

exprLaed  a'  ‘  “*  f!lght  Path  f°r  1  11"“r  ^  be 


Z(x) 


rb/2  r® 

=  /  /  f(x  "  xi<-y)  h(x  y)dx,  dy 

J-b/zJ-m  1  1 


(A. 17) 


where  b  Is  the  airplane  span  and  h(x,y)  is  the  airplane  response  along  the  flight 
path  to  a  unit  impulse  gust  impinging  on  the  wing  leading  edge  at  x  -  0  and  „  » 
Position  y  (Squat ion  (A. ,7).  as  indicated  in  Referent  Lads  L'the“  Ho,  in? 

relation  between  the  power  spectra  of  the  disturbances  and  the  response 


where 


\z  >  -  f°  I  h(Hj  ,  ar)2 )  1 2$f f  (Qi  ( n  )  d  n 

^  -CD 

'b/2  r°>  f 

lfUx,y)e'i^ix+^?  y) 


(A. 18) 


H<ni-n!>  -  r  r 

J-b/z  J -CO 


dx  dy 


3.3  Isotropic  Turbulence 

For  the  special  case  of  isotropic  turbulence 


and 

where 


Rff(£  Ay)  =  Rff/(Ax?  +  Ay2;  =  Rff(r) 


(A.  19) 
(A. 20) 


r2  =  Ax2  +  Ay2 

n2  =  +n22 

auto-correlation^funct lo^R.jf /i  xLfrom  fuTtT  “*  ‘ 

lo"ai  of 


(A.21) 


^  -  vl  AxJo^Ax)Rff(Ax)dAx 


where  J°(  )  is  the  Bessel  function  of  the  first  kind  of  order  0.  Thus  the  single 

light  runs  may  be  used  to  determine  the  spectral  characteristics  of  the  two- 
dimensional  process  f(x,y) 
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3.4  Turbulence  Components 


In  studying  the  airplane  response  to  atmospheric  turbulence,  it  is  frequently 
required  to  consider  the  effect  of  the 'several  directional  components  of  the  turbu¬ 
lence.  Thus,  f(x,y)  must  be  considered  a  vector  process  with  components  u(x,y), 
v(x,y),  and  w(x,y)  along  the  x,  y,  and  z  axes,  respectively.  The  statistical 
characteristics  defined  earlier  for  the  scalar  function  have  then  to  be  defined  for 
each  of  the  three  components  as  well  as  for  six  cross -component  terms.  These  cross¬ 
component  terms  involve  averages  of  the  form: 

Av.{u(Xj,yj)  v(x2,y2)>  =  Ruv(Ax,Ay) 


The  following  auto-covariance  tensor  covering  the  nine  functions  which  are  required 
is  of  concern: 


R(Ax,  Ay) 


Ruu  Ruv  Ruw 
^vu  ^vv  ^vw 

^wu  Rww 


For  the  case  of  isotropic  turbulence,  von  Karman  and  Howarth28  have  shown  that  these 
functions  may  all  be  expressed  in  terms  of  a  single  scalar  function  f(r)  which  depends 
only  on  the  radial  distance  r2  =  Ax2  +  Ay2.  The  relations  are 


^uu 


f(r)  -  g(r) 
r2 


(Ax)2  +  g(r) 


Ruv  =  Kyu  =  f<r)r',g(r)(A»)(Ay) 

Ryy  =  +  ,(r) 


(A.  22) 


^uw  =  “  Rvw 

^ww  =: 


~  ^wv  =  ® 


where  g(r)  and  f(r)  an.  the  covariances  between  the  velocity  components  normal  to 
and  parallel  to  the  vector  r.  The  functions  f(r)  and  g(r)  are  also  related  by  the 
following  differential  equation: 


g(r)  =  f  (r)  +  (A.  23) 

2  dr 

Thus,  for  isotropic  turbulence,  the  auto-correlation  functions  between  the  turbulence 
components  at  various  different  positions  depends  only  upon  the  single  function  f(r). 

The  relation  between  the  auto-covariances  given  by  Equation  (A. 23)  may  also  be 
used  to  derive  the  relation  between  the  power  spectra  of  the  velocity  components.  In 
terms  of  these  power  spectra,  the  relation  is  given  by 


(A. 24) 


G(0)  =  1  F(0)  -in  '^0 
2  2d 

where  G(il)  is  the  power  spectrum  of  the  lateral  gust  component  and  F(fl)  is  the  power 
spectrum  of  the  longitudinal  gust  component.  The  vertical  and  side  gu§t  components 
sensed  by  an  airplane  are  lateral  components,  while  the  gust  component  sensed  by  a 
sensitive  pitot  head  is  the  longitudinal  component.  Equation  (A. 24)  may  thus  be  used 
to  determine  one  of  these  spectra  from  measurements  of  the  other. 

Equations  (A. 23)  and  (A. 24)  are  difficult  to  apply  in  checking  on  questions  of 
isotropy  for  experimental  measurements  since  the  derivatives  of  f(r)  and  F(O)  are 
difficult  to  determine  experimentally.  However,  idealized  analytic  approximations 
to  the  measured  spectra  may  be  used  as  indicated  in  the  text  to  examine  measurements 
of  F(il)  and  G(O)  in  order  to  determine  whether  the  isotropic  condition  imposed  by 
Equation  (A. 24)  is  met  over  the  frequency  range  covered  by  the  data.  In  practice, 
it  might  be  expected  that  isotropy  will,  in  general,  apply  over  a  restricted  fre¬ 
quency  region. 
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